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We perform a combined experimental-theoretical study of the Fe-oxychalcogenides (FeOC/i) se¬ 
ries La202Fe20M2 (M=S, Se), which is the latest among the Fe-based materials with the potential 
to show unconventional high-Tc superconductivity (HTSC). A combination of incoherent Hubbard 
features in X-ray absorption (XAS) and resonant inelastic X-ray scattering (RIXS) spectra, as well 
as resitivity data, reveal that the parent FeOCh are correlation-driven insulators. To uncover micro¬ 
scopies underlying these findings, we perform local density approximation-plus-dynamical mean held 
theory (LDA-fDMFT) calculations that unravel a Mott-Kondo insulating state. Based upon good 
agreement between theory and a range of data, we propose that FeO Ch may constitute a new, ideal 
testing ground to explore HTSC arising from a strange metal proximate to a novel selective-Mott 
quantum criticality. 

PACS numbers: 74.70.Xa, 74.70.-b, 74.25.F-, 74.25.Jb 


Since the discovery [T] of high-temperature supercon¬ 
ductivity (HTSC) in iron pnictides (FePn), a fundamen¬ 
tal question has been whether the iron pnictides are 
weakly correlated metals or bad metals in close proximity 
to Mott localization [an]. These competing ideas must 
imply very different mechanisms of Fe-based HTSC. The 
Mottness view proposes that strong, orbitally-selective 
electron correlation is relevant to HTSC, while the 
competing itinerant view holds that antiferromagnetic 
(AFM) order and superconductivity (SC) originate from 
a Kohn-Luttinger-like Fermi surface (FS) nesting mecha¬ 
nism with weak electronic correlations [3] . Early insights 
suggested the importance of varying the strength HIS] 
of electron correlations using chemical substitutions to 
tune systems across band-width driven Mott localization. 
Since then, more evidence has accumulated in support of 
this proposal; FeSe is a bad-metal without AF order [B], 
while parent Fe-oxychalcogenides are electrical insulators 

[ZlIHl. 

Most FePn superconductors exhibit bad-metal normal 
states with hints of novel quantum criticality. The dis¬ 
covery of HTSC in alkaline iron selenides Ai_a;Fe 2 -ySe 2 , 
where A = K, Tl, Cs and Rb (referred to as ”245s”) inten¬ 
sified this basic debate, since SC with T^ = 30 K emerges 
from doping an AFM insulating state occuring due to Fe- 
vacancy induced band narrowing [3] . Studying more such 
materials to further develop a complete description of the 
generic electronic phase diagram of iron superconductors 
is clearly warranted. With this in mind, we investigated 
iron oxychalcogenides such as A 2 0Fe2 02M2 [A=(La, 


Na, Ba), (M=(S, Se)] [TU] because of (i) their similarities 
to LaFeOAs (fill) materials and (ii) the opportunity to 
study materials with slightly larger electron correlation 
than the pnicitides [7]. Moreover, La202Fe20M2 (M=S, 
Se) offer the attractive option to tune the strength of 
electronic correlations in one system through disorder- 
free chalcogen replacement of S by Se. The FeO Ch con¬ 
sist of stacked double-layered units La 202 and Fe 20 (S, 
Se) 2 . Fe20(S, Se )2 contains an Fe layer that is similar 
to the FeAs layer in the pnictides; in the latter, Fe^+ (d® 
configuration) is co-ordinated by O atoms but there is 
an extra Oxygen atom 0(2). The rare-earth layer also 
contains an additional 0(1) atom. Chalcogen substitu¬ 
tion results in a slightly expanded Fe-square-lattice unit 
cell of La202Fe20M2 compared to that of LaOFeAs m- 
The consequent decrease in the electronic bandwidth W 
must increase the Hubbard U, offering a disorder-free pa¬ 
rameter to subtly tune U/W across a critical value for a 
Mott insulator phase. 

In this Letter, we substantiate this reasoning by pre¬ 
senting first measurements of the unoccupied and occu¬ 
pied Fe and Oxygen density-of-states (DOS) obtained 
by soft X-ray absorption (XAS) and resonant inelastic 
X-ray scattering (RIXS). Strong Fe moment localiza¬ 
tion that can be tuned with chalcogen replacement, and, 
more crucially, the observation of incoherent Fe spec¬ 
tral weight due to Hubbard band formation marks the 
FeOCh as correlation-driven insulators. The incoherent 
spectral weight undergoes a resonant enhancement in the 
RIXS data and LHB spectral weight analysis, indicat- 
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FIG. 1: (Color online) (a) O Is XAS and (b) O Ka RIXS 
data for La202Fe20M2, (M = S, Se). RIXS intensity data 
was collected using photons with incident energies of 532 and 
534 eV. In (c) the electrical resistivity versus temperature T 
data and the simulated resistivity calculated using LDA + 
DMFT is shown. 


ing the Mott insulating character of both La 202 Fe 20 S 2 
and La 202 Fe 20 Se 2 . We discuss the implications of these 
findings for possible HTSC arising from an incoherent 
metal proximate to novel selective-Mott quantum criti¬ 
cality. 

We investigated well characterized La 202 Fe 20 M 2 
M=(S, Se) polycrystalline materials, with nominal com¬ 
positions, that were fabricated using solid state reactions 
methods described earlier [7]. High purity La 203 , Fe and 
(S, Se) powders were used as starting materials and lab¬ 
oratory X-ray powder diffraction data showed that these 
materials had minimal impurity phases; neutron powder 
diffraction data confirmed this finding. 

XAS and RIXS spectra for the FeOCh were collected 
in a polarization averaged condition. In order to re¬ 
duce the possibility of oxidation, iron oxychalcogenide 
samples were sheared under a pressure of 10“® Torr in 
a pre-chamber immediately before being placed in the 
ultra-high vacuum experimental chamber. All presen- 
tend measurements were performed in the paramagnetic 
state above 150 K. The Advanced Light Source beam¬ 
lines 7.0.1 and 8.0.1 delivered X-ray beams of 100 mi¬ 
cron spot-sizes and energy resolutions of 0.2 eV (0.2 eV) 
and 0.5 eV (0.6 eV) for oxygen(iron) X-ray absorption 
and emission, respectively. The oxygen A-edge XAS 
intensity is proportional to the unoccupied p-weighted 
DOS resulting from the Al = ±1 dipole selection rule 
for X-ray photon absorption. If the Fe-0 related bond is 
partially covalent, O 2p states hybridize with Fe 3d or¬ 
bitals and hence the pre-threshold for absorption, i.e., the 
low energy fine structure, reveals itself in ligand-to-metal 
charge transfer excitations m- XAS is sensitive to metal 




FIG. 2: (Golor online) The Fe L 2 , 3 -edge XAS profiles of 
La202Fe20M2, (M = S, Se) is shown in panel (a). The 
arrows indicate the X-ray absorption features at the photon 
energies, 707.4 and 708.8 eV, nsed to collect the Fe RIXS data 
in (b). The data are plotted with a vertical offset for clarity. 


atom coordination geometry and 3d occupation (oxida¬ 
tion state). Fig. top panel (a), shows the O K XAS 
profiles of La 202 Fe 20 M 2 (M = S, Se). Incident X-ray 
energies of 532 and 534 eV were used to collect the O Ka 
RIXS data presented in Fig. [^(a). The presence of lower 
Hubbard band (LHB) structure in La 202 Fe 20 S 2 RIXS 
spectra is consistent with resistivity data (Fig. [I](b)) 
which shows M = S to be a better insulator than M = 
Se. We note the absence of O charge transfer band, i.e., 
a high energy shoulder structure for the O 2p mainband 
peak, in the Oxygen RIXS data of La 202 Fe 20 Se 2 mate¬ 
rial: this suggests that the extent of Fe-0 hybridization 
DOS is different for M = S and M=Se compounds. 

In Fig. (a), we show the Fe T 2 , 3 -edge XAS spectra 
for both La 202 Fe 20 S 2 and La 202 Fe 20 Se 2 . A magnifica¬ 
tion (inset) of the region near the Fermi level, presented 
in Fig. (a), reveal a clear difference in the iron conduc¬ 
tion band weight of the two materials. La 202 Fe 20 Se 2 
has enhanced spectral weight at pre-edge just above Ep, 
and this must ultimately be related to details of the cor¬ 
related electronic structure changes upon replacement of 
S by Se. Fe XAS data is consistent with both RIXS 
and bulk measurements that show M = Se to be less 
strongly correlated than M = S. Resonant inelastic X-ray 
scattering (RIXS) intensity data, see Fig. [^(b), were col¬ 
lected using incident X-ray energies tuned to near- and 
off-resonance Fe X-ray absorption features at 707.4 eV 
and 708.8 eV, respectively. RIXS [H] intensity results 
from a second-order (i.e., two-step) photon-in, photon- 
out process that can be can be described by the Kramers- 
Heisenberg (KH) relationship [12]. The KH expression 
requires the enhancement of the scattering intensity when 
incident photon energy matches X-ray absorption edge 
energies. The RIXS data for M = S shows correlation 
induced broad features in the high energy spectral re- 
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FIG. 3: (Color online) Orbital-resolved LDA density-of-states 
(DOS) for the Fe d orbitals of La202Fe20Se2 as well as 
LDA-I-DMFT results for different values of U (with U' = 
U — 2Jh) and fixed Jh ~ 0.7 eV. Notice the narrow bands in 
the LDA DOS. Compared to the LDA results, large spectral 
weight transfer along with electronic localization is visible in 
the LDA-I-DMFT spectral functions for 17 > 4 eV. 

gion that form the LHB. In the RIXS data, for both 
M = S, Se, there is a high intensity Fe 3d metal va¬ 
lence band peak. Near 692 eV, a broad shoulder in the 
high energy spectral region is resonantly enhanced. The 
spectra clearly reveal incoherent electronic excitations 
involving 3d Fe electrons. These are correlation satel¬ 
lites formed by multi-particle excitations that involve lo¬ 
calized states, and clearly reveal the stronger electronic 
correlations that generate a correlated insulator in the 
FeO Ch. 

Taken together, these results mark the FeO Ch as corre¬ 
lated, bandwidth-control driven insulators m- In multi¬ 
orbital FeO Ch with strong inter-orbital charge transfer, 
correlations and crystal field effects, proper characteriza¬ 
tion of the Mott insulating state(s) is an important, non¬ 
trivial issue. To shed light on microscopies underlying the 
above findings, we performed LDA-l-DMFT calculations 
for the FeO Ch closely following earlier work PE]. LDA 
calculations were performed within the linear muffin- 
tin orbitals (LMTO) scheme. The one-electron Hamil¬ 
tonian reads Hq = J2k,a,a with a = 

xy, xz, yz, x^ — y‘^,3z‘^ — label (diagonalized in orbital 
basis) five Fe 3d bands. We retain only the d states, 
since the non-d orbital DOS have negligible or no weight 
at Ep. An important aspect of our study is readily vis¬ 
ible in Fig. which shows a sizable reduction of the 
LDA bandwidth relative to that of FeSe (which is it¬ 
self an incoherent bad-metal m), induced by hybridiza¬ 
tion with O states and distorting effects of the La-layers. 
Another interesting aspect, in contrast to Fe-arsenides, 
is that the xz — yz orbital degeneracy in the tetrago¬ 
nal structure is explicitly removed at the outset. Thus, 
in FeO Ch, an electronic nematic (EN) instability linked 
to ferro-orbital order and its potentially related critical¬ 


FIG. 4: (Color online) Orbital-resolved imaginary parts of the 
DMFT self-energies for the Fe d orbitals of La202Fe20Se2 for 
different values of U and Jh ~ 0.7 eV. 


ity is not relevant. Finally, clear lower-dimensional band 
structural features are visible in Px,^^(w) (quasi-lD) and 
pTda(^) (2D van-Hove singular) near the Fermi energy, 
Ep{= 0). Enhancement of the effective U/W ratio in 
such an intrinsically anisotropic setting should now nat¬ 
urally favor the Mott insulator state in FeO Ch. To sub¬ 
stantiate this, DMFT calculations were performed us¬ 
ing the multi-orbital iterated perturbation theory (MO- 
IPT) as an impurity solver. Though not exact, MO-IPT 
has a proven record of providing reliable results for a 
number of correlated systems m, and is a very effi¬ 
cient and fast solver for arbitrary T and band-hlling. 
The full Hamiltonian is H = Hq + Hi with Hi = 

with 

U' = {U—2Jh) being the inter-orbital Coulomb term and 
Jh the local Hund’s coupling. Though U ~ 4.0 — 5.0 eV 
and Jh = 0.7 eV are expected, we vary 3.0 < D < 5.0 to 
get more insight. 

In Fig. P we show the total many-body local spectral 
function (DOS) for the M = Se case. A clear transition 
from a bad-metal state at U = 3.0 eV to a correlated 
insulator for U = 4.0, 5.0 eV occurs. Looking closely at 
the self-energies in Fig. however, a remarkable aspect 

stands out: ImSa(w = Ep) vanishes in the region of 
the Mott gap, instead of having a pole, as would occur in 
a true Mott insulator. This aspect is more reminiscent of 
a correlated Kondo insulator, where the gap arises due to 
combined effects of electronic correlations and sizable in¬ 
terband hybridization. While this novel finding is some¬ 
what unexpected, proximity to a true Mott state is seen 
by observing that the pole structures in ImEa(a;) with 
a = xz,yz,3z^ — lie only slightly above Ep. Thus, 
because the system is correlated, sizable spectral redis¬ 
tribution in response to additional small perturbations 
(e.g, longer-range Coulomb interactions in an insulator) 
can drive it into a true Mott insulating phase. More 
importantly, electron doping will also generically cause 
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FIG. 5: (Color online) Total LDA+DMFT XAS spectrum of 
La20Fe202M2 (M =S, Se), showing the changes in the up¬ 
per Hubbard band upon chalcogen substitution. Inset shows 
the corresponding total LDA-I-DMFT spectral function for 
U =5.0 eV and Jh =0.7 eV. 


large-scale orbital-selective (OS) spectral weight trans¬ 
fer, leading to orbital-selective Mott transition (OSMT), 
wherein a subset of orbital states remain (Mott) gapped 
in an incoherent metallic state [5]. Such states can ex¬ 
hibit novel quantum criticality associated with an end¬ 
point of the OSMT. Whether HTSC near such novel crit¬ 
icality occurs in suitably doped FeO Ch is an enticing and 
open issue. 

Remarkably, a comparison to Fe-edge XAS and dc 
resistivity data reveals good agreement, providing sup¬ 
port for this novel finding. In Fig. we show our 
LDA-I-DMFT XAS spectrum. We have incorporated 
the relevant effects of core-hole scattering for XAS line- 
shapes by consistently adapting the procedure of Pardini 
et al. m within our DMFT scheme. Several important 
aspects stand out: (i) there is good semi-quantitative 
accord with with the Fe-edge XAS data up to 3.0 eV 
above Ep. More importantly, details of the differences 
in spectral weight between M =S and M =Se cases is 
also in very good accord with data, (ii) enhanced spec¬ 
tral weight at pre-edge just above Ep, just as seen, while 
smearing of lower energy structures (0.0 < a; < 2.0 eV) 
and higher energy of the 3.0 eV peak for M =S relative to 
M =Se is in good accord, both with XAS data, as well 
as with resistivity data which show M =S to be more 
correlated than M =Se. Indeed, LDA-I-DMFT resistiv¬ 
ities (using the Kubo formalism and neglecting vertex 
corrections to the conductivity, which are negligible m) 
for M =S, Se show very good accord with experiment 
(where was measured above 150 K. This is a 

strong internal self-consistency check, and puts our novel 
proposal on stronger ground. 

Armed with these positive features, we now discuss 


their implications for possible HTSC in suitably doped 
or pressurized systems. Both are expected to tune the 
Fe-0 hybridization: since we find a correlated Mott- 
Kondo insulator, this must non-trivially reconstruct the 
Fe-d states. In La202Fe20Se2, suppression of oxygen 
Hubbard-band spectral weight presages the onset of bad- 
metallicity even as the localization of the Fe id states is 
preserved. This may provide clues to metallize YeOCh 
by appropriately manipulating the Oxygen and chalcogen 
states which hybridize with Fe-d states. This is an ex¬ 
citing perspective because, if HTSC results, (as in 245s), 
it would reinforce the fundamental challenge to the itin¬ 
erancy perspective in pnictide HTSC research. Specifi¬ 
cally, since strong correlations have now obliterated the 
LDA FS as above, large-scale orbital-selective redistribu¬ 
tion of spectral weight upon carrier doping is generically 
expected to lead to OSMT and FS involving only a sub¬ 
set of original d-orbitals. The LDA FS and its nesting 
features can now no longer serve as a guide for rational¬ 
izing instabilities to novel order(s). Instead, strong elec¬ 
tronic scattering between the partially (Mott) localized 
and itinerant subsets of the many-body spectral func¬ 
tions generically extinguishes the Landau quasiparticle 
(QP) pole at Ep. In absence of a QP description for 
the normal state, the very tenability of the FS nesting 
scenario is called into doubt. 

Thus, our work provides an impetus to consider closer 
similarities between Fe-based systems and the cuprates. 

In similarity with 245 materials, we suggest that 
FeO Ch are ideal candidates for testing this proposal. 
The development of incoherent scattering observed in 
our XAS and RIXS spectra arises from Mott localized 
states. If these can be appropriately doped, it is likely 
that strong orbitally-selective scattering will nearly ex¬ 
tinguish the Landau-Fermi liquid quasiparticles, leading 
to the emergence of an incoherent pseudogapped metal 
reminiscent of underdoped cuprates [H]. Whether novel 
quantum criticality associated with a selective-Mott tran¬ 
sition in conjunction with the development of HTSC in 
FeO Ch will be seen in the future is of great interest. Our 
findings provide a compelling motivation to produce and 
study doped or pressurized iron oxychalcogenides. 

After the recent journal submission of this paper, we 
learned of work by G. Giovannetti et. al. on the iron 
oxychalcogenides which is complementary to ours and 
gives comparable results. 
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